en 

o 

(N 

Oh 
< 

0^ 



I 



> 

(N 

in 
in 

o 



An Infinite Set of Ward Identities for Adiabatic Modes in Cosmology 

Kurt Hinterbichler", Lam Hui^ and Justin Khoury'^ 



" Perimeter Institute for Theoretical Physics, 
31 Caroline St. N, Waterloo, Ontario, Canada, N2L 2Y5 



^ Physics Department and Institute for Strings, Cosmology and Astroparticle Physics, 
Columbia University, New York, NY 10027, USA 



Center for Particle Cosmology, Department of Physics & Astronomy, University of Pennsylvania, 

209 South 33rd Street, Philadelphia, PA 19104 



Abstract 



^ ^ We show that the correlation functions of any single- field cosmological model with constant growing- 

s' modes are constrained by an infinite number of novel consistency relations, which relate N + 1-point 

^ correlation functions with a soft-momentum scalar or tensor mode to a symmetry transformation on A'^- 

point correlation functions of hard-momentum modes. We derive these consistency relations from Ward 
identities for an infinite tower of non-linearly realized global symmetries governing scalar and tensor 
perturbations. These symmetries can be labeled by an integer n. At each order n, the consistency 
relations constrain — completely for n = 0, 1, and partially for n > 2 — the q" behavior of the soft 
limits. The identities at n = recover Maldacena's original consistency relations for a soft scalar and 
tensor mode, n = 1 gives the recently-discovered conformal consistency relations, and the identities for 
n > 2 are new. As a check, we verify directly that the n = 2 identity is satisfied by known correlation 
functions in slow-roll inflation. 



1 Introduction 

The consistency relation is one of the most powerful probes of early universe physics. It states that the 
3-point function of the curvature perturbation (^ [1, 2], in the squeezed or local limit where one of the 
modes becomes soft, is determined by the scale transformation of the 2-point function^ [3~5]: 



1 , - 9 



li- ^C^^.U)' = -^1 • ^^^^aC.^.)' ■ (1-1) 



Combined with the observation that the 2-point function is scale invariant to within a few percent, this 
tells us that the 3-point function in this limit should be vanishingly small. The consistency relation 
holds under very general conditions [5]: any early universe scenario involving a single scalar degree of 
freedom (or single 'clock') whose perturbations grow to a constant at late times, must satisfy (1.1). 
This encompasses all single-field inflationary models [6-9], including hybrid models [10], as well as some 
non-inflationary scenarios [11-21]. Conversely, to observe a significant 3-point function in this limit 
requires violating one of the assumptions above, for instance through additional scalar fields [22-27] or 
an unstable background [28, 29, 32, 33]^. The consistency relation (1.1) has been established through 
explicit calculations [3, 5] and can be understood using "background- wave" arguments'^ [3, 4]. It holds 
whether the short-wavelength modes are inside or outside the horizon [37]. 

In this paper we will show that inflationary correlation functions — and more generally those of any 
single-field cosmology with constant growing mode solutions for scalars and tensors — are constrained 
by an infinite number of consistency relations. Similarly to (1.1), they relate the soft limit g — >• of a 
scalar or tensor mode in an A^ -|- 1-point correlation function to a symmetry transformation on an A^-point 
function. We will show how these arise as Ward identities for an infinite number of non-linearly realized 
symmetries, akin to the soft-pion theorems of chiral perturbation theory [38, 39]. At each order n, the 
consistency relations constrain — completely for n = 0, 1, and partially for n >2 — the g" behavior of 
the soft limits of correlation functions. 

As the simplest example, we will recover the original consistency relation (1.1) from the Ward 
identity for spontaneously broken spatial dilations. This was already pointed out in [3] and recently shown 
in [40]. The operator approach followed here in translating the Ward identity to a consistency relation has 
some overlap with [40]. The approach establishes the non-perturbative nature of the consistency relations. 
In a parallel investigation, Goldberger, Hui and Nicolis [41] have used the path integral approach to 
derive (1.1) from a Ward identity. See [42] for a recent derivation of (1.1) using holographic arguments. 
As another example, we will recover the tensor consistency relation [3] 

^IW)^^^^"^^'^"^'^^' = -r^^o('/)^l^(Q,C,^.)'. (1.2) 



""^Here P,^ denotes the power spectrum, and (. . .)' are correlators without the momentum-conserving 5-function. 

■^One exception, for instance, is the Solid Inflation model proposed recently [34]. 

^The background- wave argument basically goes as follows. The long ((f — >■ 0) mode freezes out much earlier than the short 
(^1,2) modes, and hence acts as a classical background for the generation of these modes. Since the growing mode solution 
for 1^ is a constant, the background 3-metric (ignoring tensors) experienced by the short modes is hij = e^''^ a? {t)5ij , where 
("l is approximately constant. This constant background is an adiabatic mode [35, 36], which can be removed by a local 
rescaling of the short modes. 



where ef^ is the polarization tensor, from the Ward identity for spontaneously broken anisotropic spatial 
rescaling. 

Recently, it has been pointed out [43, 44] that scalar perturbations in spatially flat single-field 
cosmology are governed by the symmetry breaking pattern 

SO{A, 1) — )• spatial rotations + translations , (1-3) 

with C playing the role of the Goldstone boson (or dilaton) for the 5*0(4, 1) conformal symmetries on 
M?. The origin of conformal symmetry is most easily seen in comoving gauge, where the scalar field 
is unperturbed, and the spatial metric (ignoring tensors) is hij = a^{t)e'^^^^'^'5ij. Since this 3-metric is 
conformally flat, the 10 conformal transformations on K,^ preserve this form and hence are symmetries of 
the scalar sector."* The dilation and special conformal transformations (SCTs) are non-linearly realized on 
C, whereas spatial translations and rotations are linearly realized. The conformal symmetries are restored 
(and linearly realized on the fields) in the limit that the background becomes exact de Sitter space [45-52]. 

As another special case of our general Ward identities, we will show that the conformal consistency 
relation, recently derived using background- wave arguments [43], is a consequence of the Ward identity for 
non-linearly realized SCTs. The conformal consistency relation relates the order q behavior of correlation 
functions with a soft-<^ mode to an SCT acting on the correlation functions of the hard modes. Since 
SCTs result in a departure from the transverse, traceless conditions for tensor modes, they are strictly 
speaking a good symmetry of the scalar sector only [43, 44]. However, we will show how the SCTs can 
be corrected, order by order in the tensors, to become full-fledged symmetries of the theory of scalar 
and tensor perturbations. These result in corrections to the conformal consistency relation to all orders 
in tensors. (In contrast, the spatial dilation symmetry is exact, including tensors, and (1.1) receives no 
correction.) There is also a tensor consistency relation constraining the order q behavior of the correlation 
functions with a soft tensor mode [43], which we will also recover from our Ward identities. 

More generally, we will show in Sec. 2 that the theory of scalar and tensor perturbations in comoving 
gauge, 

/i,, = a2(t)e2«-^'*) fe^(-^'*)) , (1.4) 

where^ 7*^ = 0, 9i7* • = 0, is constrained by infinitely-many non-linearly realized symmetries. These are 
residual diffeomorphisms, which do not fall off at infinity and which leave the form of the 3-metric (1.4) 
intact. The symmetries are defined perturbatively in the tensors, and the corresponding field transfor- 
mations include terms to all orders in 7. No quasi-de Sitter or slow-roll approximation will be assumed 
— the symmetries uncovered here hold on any spatially-flat cosmological background. 

Since the residual diffeomorphisms diverge at infinity, they map field configurations which fall off 
at infinity into those which do not. Nevertheless, in Sec. 2.3 we will generalize Weinberg's argument [36] 
to show that certain linear combinations of these transformations can be smoothly extended to physical 

^This does not contradict the usual notion that "comoving gauge completely fixes the gauge," since this statement assumes 
that gauge transformations fall off at infinity — conformal transformations clearly do not. As such, conformal transformations 
are residual diffeomorphisms mapping field configurations which fall off at infinity into those which do not. Nevertheless, as 
shown in [43, 44] and reviewed in Sec. 2.3, they can be extended to transformations which do fall off at infinity, and hence 
generate new physical solutions or adiabatic modes [35, 36]. 

^Unless otherwise stated, spatial indices are everywhere raised/lowered using 5ij. 



configurations which fall off at infinity. In other words, they correspond to adiabatic modes. Such 
transformations can be thought of as the g — )• limit of transformations which do fall off at infinity, and 
therefore generate new physical solutions. 

In Sec. 3 we will show how the Ward identities associated with the non-linear ly symmetries lead to 
consistency relations for the soft limits of correlation functions. The derivation, while somewhat technical, 
is very general and allows us to derive at once all consistency relations. (In contrast, the background- 
wave method, while straightforward for the dilation consistency relation [4], is already considerably more 
intricate for the conformal consistency relation [43].) Another upshot of the operator approach is that 
the regime of applicability of the consistency relations is sharply defined — Ward identities are non- 
perturbative statements, and in particular do not rely on semi-classical approximations. 

The consistency relations we obtain are of the following schematic form (with n > 0) 

where O represents a product of N scalar (" and tensor 7 perturbations with momenta labeled hy ki, ..., k^; 
O represents the same set of fields with one of the fields (C or 7) replaced by another 7; ( )^ represents 
the connected momentum space correlation functions with the overall delta function removed. On the 
left hand side, we have N + 1-point functions with a soft (small q) C ^"^ 1 l^g (some relations have only 
soft C,, some only soft 7, and some have both); on the right hand side, we have A^-point functions with 
the soft leg removed. As such, these consistency relations resemble the well-known soft-pion theorems, 
and indeed they follow from Ward identities applied to nonlinear ly realized symmetries just like the soft- 
pion theorems do. There is an infinite set of consistency relations, labeled by the integer n > 0, each 
controlling the q^ behavior of the relevant A^ + 1-point function. The other momenta fci, ..., /ctv are hard 
compared to q, but there is no assumption on whether they are inside or outside the horizon [37]. There 
are 3 independent relations for n = (one involving a soft scalar and two involving a soft tensor), 7 
relations for n = 1 (three involving a soft scalar and four involving a soft tensor), and 6 for each n > 2 
(four involving a soft tensor and two involving mixtures of soft scalar and tensor). The n = and n = 1 
relations are known. The n > 2 relations are new. 

In Sec. 4 we will show that the n = identities reproduce the dilation and anisotropic scaling 
consistency relations, of which (1.1) and (1.2) are the simplest renditions. Moreover, we will recover from 
the n = 1 identity the conformal consistency relation and the linear-gradient tensor analogue relation. 
While the n = 0, 1 consistency relations completely fix the q^ and q behavior of the soft limits, the 
n > 2 identities only partially fix the (7" behavior of the correlation functions. Of the 6 algebraically- 
independent consistency relations at n > 2,4 involve soft-7 correlation functions only, while the remaining 
2 involve a linear combination of soft-C and soft-7 correlators. (This is a key difference compared to the 
lower-order identities — there are no "pure scalar" consistency relations for n > 2.) As an example 
of a novel consistency relation, we study in Sec. 5 the n = 2 tensor Ward identity relating the q^ 
behavior of (7ij(Q)Cg C^ ) in the limit g* — )• to derivatives of the scalar 2-point function ((^g C^ )■ ^^ 
a check of our result, we verify that this new consistency relation is satisfied by the 3-point function 
computed by Maldacena in slow-roll infiation [3]. We conclude in Sec. 6 with a discussion of possible future 
investigations. Many of the technical details of our derivation are described in a series of Appendices. 



The general Ward identities derived here constitute the complete list of consistency relations that 
any single-field cosmological scenario with constant growing modes must satisfy. The Planck data, which 
has found no evidence of non-Gaussianity in the squeezed limit [53], is so far consistent with the lowest- 
order consistency relation (1.1). The higher-order relations involving scalars and tensors, while unlikely 
to be tested in the near future, offer the complete checklist with which current and future-generation 
observations will put single-field inflation to test. 

2 Symmetries of Cosmological Perturbations 

Consider perturbations on a spatially-fiat, Friedmann- Robertson- Walker (FRW) background driven by a 
scalar field (f) = (j){t). Following [3], we work in comoving (or uniform-density) gauge^, where the scalar 
field is unperturbed, 

ck = m. /i^, = a2(t)e2';(-^'*)(e^(-^'*))^^, f, = , 9^7^=0. (2.1) 

In this gauge, scalar inhomogeneities are captured by Q, while tensor modes are encoded in 7jj. 

If all the fields and gauge parameters are assumed to die off sufficiently fast at spatial infinity, then 
this gauge choice completely fixes the gauge, leaving no residual symmetries. Nevertheless, there is still 
room for diffeomorphisms whose gauge parameters do not die off at infinity. In the absence of tensors, 
for instance, it has been argued recently that conformal transformations on spatial slices are residual 
symmetries of the scalar sector [43, 44]. Our goal in this Section is to generalize these results to include 
tensor perturbations to all orders. Although we will eventually be interested in deriving Ward identities for 
inflationary correlation functions, we should stress that no slow-roll or quasi-de Sitter approximations will 
be necessary to identify these residual symmetries. They hold on any (spatially-flat) FRW background. 

Consider a (possibly time-dependent) spatial diffeomorphism, ^*(x, t). Since this diffeomorphism is 
purely spatial, it does not disturb the gauge choice (p = (p{t). If we can identify fleld transformations 5^ 
and 5jij such that 

^ (e'^ i^'h) = A (^'^ i^'h) ' (2-2) 

then these transformations will act like diffeomorphisms which leave the gauge fixed action invariant, and 
hence constitute a symmetry. The Lie derivative C^{gij) is defined as C^d^gij +^^digkj +S,^djgik- We will 
solve (2.2) order by order in powers of 7, expanding the field variations and diffeomorphism parameter 
as follows: 

At each order we will impose that 5jij be transverse and traceless. 



i. = «l'°'+?f ■' + ... (2,3) 



®The form of hij in this gauge is closely related to the conformal decomposition introduced by York [54] and Lichnerow- 
icz [55]. 



2.1 Zeroth-order in tensors 

We first focus on the zeroth-order diffeomorphisms, ^'^ , leaving the derivation of the first-order correction 
e-^'^ to Sec. 2.2 below. 



At zeroth-order in 7, expanding (2.2) gives 

2S&'% + S,f^ = 2cf^d\5., + d^'^ + a,ej^") . (2.4) 

Since 6^ij must be traceless, taking the trace allows us to solve for 5C ' , 



6C^-'^ = ld^^f^+cf^dX. (2.5) 



Plugging back into (2.4), we can solve for 6^1- , 



,,(/) = 8.5<^»)+a,£r)-?8-{r%«. (2,6) 

Taking the divergence of this expression yields an equation for Q"' , 

V'('f"^ + l^^^^^f^=0■ (2.7) 

Note that (2.7) admits no non-trivial solutions which vanish at infinity. We see this by taking its diver- 
gence, V'^{d'^Q"' ) = 0, which for fields that vanish at infinity implies d'^Q"' = 0. Plugging this back 
into (2.7) gives V^^^^''^ = 0, which then implies ^'^ = 0. Hence Q'^ cannot vanish at infinity. As a 
corollary, even though we allow for q"' to depend on C, the argument above implies that it is in fact 
independent of C- To see this, imagine expanding Q in powers of (, where each term in this expansion 
must separately satisfy (2.7). Since higher-order terms are proportional to powers of C, and thus vanish 
at infinity (since C does), they must be trivial. 

Any (possibly time-dependent) q"' satisfying (2.7) preserves comoving gauge, to zeroth order in 
tensors, with field transformations (2.5) and (2.6). The condition (2.7) is the divergence of the conformal 
Killing equation on R'^, di^j + dj^i = ^dk^ Sij, so any of the transformations of the conformal group 
on R'^ (with possibly time dependent parameters) will satisfy it. These include, in addition to (time- 
dependent) spatial rotations and translations, the spatial dilation and 3 special conformal transformations 
(SCTs) [43, 44] 

^dilation ^ ^(^)^. 

^fCT ^ 2h> {t)xjX^-xHi{t). (2.8) 

2.2 Higher order in tensors 

It is worth emphasizing that while SCTs were initially believed to be good symmetries of the scalar 
sector only [43, 44], since they result in a departure from the transverse, traceless conditions for tensor 
modes, here we show that ^f can be corrected order by order in 7 to preserve the transverse, traceless 



conditions, and thus is promoted to a full-fledged gauge-preserving transformation. In contrast, we also 
show that spatial dilation is exact to all orders in tensors. 

At first order in 7, (2.2) gives 

2<5C(^^)<5., + 25C(^'\., + 6^1]'^ = 2^'^'^^'C6^, + 2^'^'^^''C7^J + d^^'^ + d,C^'^ + ^^(.0)7., • (2.9) 

The tensor variation 57^.^ must be transverse and traceless, so as before, taking the trace allows us to 
solve for the scalar variation, 

S^h') = ^^')dX + \nt'^^ + \d.ifh'' . (2.10) 

Plugging back into (2.9), we can solve for the tensor variation, 

-^8i4"°'t"'5« - Is'tt^'hi, + £jh«,7.i . (2.11) 

Taking a divergence, we find the following equation for Q , 

^Hf' + \d. {d^^^'^) = d^ [id,ct'y% + id'^t'y. - ^,(.0)7..) . (2.12) 

The right-hand side contains known zeroth-order quantities only and thus sources Q"' . We are instructed 
to solve this equation subject to the boundary condition that Q^ vanish at infinity. (Since Q^ is first- 
order in 7, which itself goes to zero at infinity, there is no room for a homogeneous solution that would 
not be proportional to 7.) By taking a divergence of (2.12) it is easy to see that the solution is unique, 
albeit spatially non-local. Explicitly, the solution for Q with these boundary conditions is 

^P'' = ~ (^/ - ~) (^,(.0)7H - \d^ifh-'5^' - Id-i^^y^') . (2.13) 

There is no obstruction to extending this to any order in 7. At m-th order in the tensors, (2.2) 
reads 

25C^^"'h,j + 5jif^ = 2^l?"^'^d''C6ij + 9iCJ^'"^ + dj^i^"^^ + (known lower order pieces) . (2.14) 

As before, taking the trace allows us to solve for 6Q'^"^\ and then plugging back we can solve for d^^^'J . 

f'7"^") 

Taking the divergence then yields an equation for Q which reads 

V Q'^ + -di idj^^^"' ' ] = (known lower order pieces) . (2.15) 

The same line of argument as before leads us to conclude that (2.15) has a unique solution. 

(7") dil 

Note that the dilation transformation, corresponding to Q ' = Xxi, is exceptional in that it 

does not get corrected at any order in the tensors. Indeed, from (2.13) it is straightforward to show that 
^j- ' ■ = in this case. This is true to all orders: ^f^^' = Xxi is an exact symmetry, with 

^dilation^(^) = A (1 + x'd.ax)) , 
^dilation^.^.(^) = Xx'^dkJijix)- (2.16) 



2.3 Adiabatic modes 

By the argument above, the residual transformations S,i have a field-independent part which does not fall 
off at spatial infinity, and therefore map field configurations which fall off at infinity into those which 
do not. Nevertheless, by generalizing Weinberg's adiabatic mode argument [36], we will see — at least 
to linear order in perturbations — that a subset of these transformations can be extended to physical 
configurations with suitable fall-off behavior at infinity. In other words, this subset of transformations 
can be thought of as the g — >• limit of transformations which do fall off at infinity, and therefore generate 
new physical solutions. 

To be extendible to a physical mode, a configuration must satisfy the equations of motion away 
from zero momenta, that is, it cannot "accidentally" solve the equations simply because it is being hit by 
spatial derivatives. The only equations for which this may happen are the constraint equations of General 
Relativity, since the evolution equations have terms with second time derivatives and hence no spatial 
derivatives (since all the equations are second order). At linear order, the momentum and Hamiltonian 
constraints are 

2di (^HNi - C) - ^V^iV, + ^di {djN^) = ; 

di {dX + HN') + ^Ni + Sa^H (hNi - c) = . (2.17) 

Cg 

Here, A'^ and Ni represent the lapse and shift of the metric, Ni is the perturbation of A^, H is the Hubble 
parameter, and Cg is the sound speed. For field configurations with suitable fall-off behavior at spatial 
infinity, we can solve these equations for A^i and A^* by inverting V^ following [3]: 

Meanwhile, to linear order we only need the non-linear part of the transformation laws, 

SC = \die, SN' = e, SN,=0. (2.19) 

Since this represents a diffeomorphism from the unperturbed solution, and the equations of General Rel- 
ativity are diffeomorphism invariant, these transformations automatically preserve the constraints (2.17). 
To be extendible to a physical field configuration, however, it must also be consistent with the solu- 
tion (2.18). To start with, since 6Ni = 0, the first of (2.18) implies 

diC = . (2.20) 

Meanwhile, the solution for A^* requires 

r = -^5.5,e^'- (2.21) 

Any solution to (2.20), (2.21) and the gauge-preserving condition (2.7) represents a gauge transformation 
of the unperturbed solution which can be extended to a physical field configuration, i.e., an adiabatic 
mode. 



V^^i + ^did^Cj = . (2.26) 



Let us decompose the difFeomorphism as 

e = e+c^, (2.22) 

where d^^f = 0. While not unique (since ^* does not faU off at infinity), such decomposition can always 
be done. Prom (2.20), we conclude that 

diC = . (2.23) 

Any time-dependent contribution to ^' must be divergence- free, and hence can be absorbed into a redef- 
inition of ^?Y- We can therefore assume ^* is time-independent without loss of generality. Equation (2.21) 
then reduces to 

iT = -^didje, (2.24) 

with solution 

Finally, to satisfy (2.7) at all times, ^* must itself be a solution, i.e., 

1 
3 

The physically allowed diffeomorphisms are therefore given by 

where ^j is any time- independent, gauge-preserving transformation, that is, any time-independent trans- 
formation satisfying (2.26). For example, the physical extension of a time-independent SCT ^ = 
2WxjXi — x^hi, which as mentioned earlier satisfies (2.26), is 

ii = 2Wx,x,-x\-2bi -— -, (2.28) 

where the correction is recognized as a time-dependent translation, in agreement with [43, 44]. A time- 
independent spatial dilation, ^^iiation _ ^^j^ ^^ ^-^q other hand, does not get corrected and therefore 
corresponds to an adiabatic mode [36]. 

More generally, note that Q always transforms linearly under the time-dependent correction ^^, 
since 5j^^ = by definition. Specifically, the field transformations (to zeroth order in tensors) are 

1 - - ^ - - /"* dt' 

K=^die + iidX; 6jij = Hj+^^di^j + dj^i) J j^y (2.29) 

where, according to (2.6), 6jij = di^j + dj^i — 2d^^k^ij/3. 

If ^* is itself transverse such that the two terms of (2.26) vanish separately, then the physical 
symmetry receives no time-dependent correction: 

r = el, (2.30) 



with Sj^^ = V^^^ = 0. We will refer to such symmetries as "tensor symmetries", in the sense that they 
effect nonlinear shifts in the tensor but not the scalar perturbations. 

The constraints that we derive in this Section, that our residual diffeomorphism in fact solves the 
constraints at finite momenta, are necessary only for the 0(7*^) part of the diffeomorphism (we adopt the 
implicit notation that ^i without superscript "7''" is the 0(7^) part). For the 0(7^) and higher parts of 
the diffeomorphism, since 7 falls off at infinity, there is no need to check that the constraints are satisfied, 
and we simply extend the 0(7^) part of the physical transformation to higher order in the tensors using 
the procedure of Sec. 2.2. 

2.4 Taylor expansion 

As shown in (2.27), the physical symmetries are expressed in terms of time-independent diffeomorphisms 
^* satisfying the gauge preserving condition (2.26). In general, this can be solved in a power series: 



00 00 ^ 



(2.31) 



n=0 



n=0 



where the array Mj^^...^^ is constant and symmetric in its last n + 1 indices.^ Equation (2.7) translates 
to the condition: 

Mim2...e,. = -^Muu,..!^ (n > 1) . (2.32) 



Each diffeomorphism ^^ in the series generates its own field transformations (2.5) and (2.6). To derive 
the Ward identities consistently at lowest non-trivial orders in the tensors, we will need to work to linear 
order in 7 for the field transformations. Using (2.13), the transformations in momentum space for each 
n are given by® 



SC^^^k) 



+ 



ht'Xk) 



nj 



+ 



3n! 



-M 



gn 



ii£l...£„ 



-M,- 



n\ 



-^^^Hio...in 



dki^ ■ ■ -dk, 



{{2^f5\k) 



+ 



]Cik) 



n\ 



-M«„..,,„T««^^-(fc) 



dkg^ ■ ■ ■ dki^ n + 1 dk^^ ■ ■ ■ dki^ j 

QU 



dki-^ ■ ■ • dk£ 



-lij{k) + ... 



n\ 



M,-,.,...,„ + M,,,„..,„ - -6ijMeei,...i^ J ^^^^ ^^^^ 



gn 



{(.^^fsHk) 



^M«„...,„ ( 5«o ^" 



+ 



k' 



Qn+l 



dke^ ■ ■ ■ dki^ n 



-i)" - (9" 

— M,,,,.„,„rv,(fc) 



+ 1 dk(,Q ■ ■ ■ dki 
l^'{k) + ... 



\lij{k) 



(2.33) 



^We will ignore the constant term in the Taylor series expansion, ^i = Mi, since this represents an unbroken (time- 
independent) spatial translation. 

*See Appendix A for Fourier transform conventions. Our symmetrization convention is T^ab) = \ [Tab +31,a)- Note that 
the transformations on ( and 7 expressed here are affected by the time- independent diffeomorphism ^\ not its time-dependent 
deformation (into an adiabatic mode). 



9 



where 

1 1 , . 5 - - 

'^rsijM(k) = 2 \5s(^i - krk{i\ 5j)(fc<J£)r ~ \^ij ~ hkj\ Sj.{k^i)s - ^ '^J(fe^%'^rs , (2.34) 

with k^ = k^/k. The eUipses in (2.33) indicate higher-order corrections in the fields. Specifically, in 6^ 
the ellipses include terms of order (^7, 7^, etc., while in ^7 they include terms ~ 7^ and higher-order in 7. 

One further "adiabatic" condition must be imposed for these transformations to correspond to 
physically-realized symmetries. As before, we want to consider only those configurations which can be 
smoothly extended to (i.e., thought of as the long- wavelength limit of) a physical mode with suitable 
fall-off behavior at spatial infinity. The additional condition arises by demanding that the non-linear shift 
in 7j£g in (2.33) remains transverse when extended to a physical mode. To see the constraint, we imagine 
smoothing out the momentum profile around q = 0. To ensure that transversality is preserved in Fourier 
space at finite momentum, St^lUoi'l) — 0' ^^ ^^^ ^^^ ^ihr-^n coefficients become g-dependent such that^ 

M,^^i,„jSq) + Mi,,^,,„tM - \5ri,Mm,..iMj = ■ (2.35) 

To determine the space of possible M's, we proceed as follows: we first determine the M's for some 
reference direction in momentum space, say g = 2 = (0, 0, 1). These are all arrays which satisfy (2.35) for 
q = z as well as (2.32), and which are symmetric in the last n + 1 indices. We then get the remaining q 
dependence by applying to all the indices of M some standard rotation which takes z into q. 

To illustrate the construction, we discuss the lowest-order transformations in some detail: 

• n = 0: In this simplest case, the symmetry transformation is linear, Q^~ = Mj£(,x '^. There are no 
constraints coming from symmetry in the last indices or from (2.32). A general array Mj£g can be 
decomposed into a trace part, a symmetric traceless part and an anti-symmetric part: 

Mj^P = A5i£(, + Siig + ujiig , (2.36) 

where S^^^ = S^^^i, cong = —uji^i, and Su = 0. The anti-symmetric part, parametrized by uja^, 
corresponds to spatial rotations, which are linearly realized and hence will not concern us. The trace 
part, parametrized by A, corresponds to a dilation, Mj^jJ*^*'""^ = X6iig. Finally, the symmetric traceless 
part, parametrized by Sn^, describes a volume-preserving, anisotropic rescaling of coordinates, 
M^^^^^ = Siig, under which ( transforms linearly and 'jUg shifts by a constant. The constraint (2.35) 
tells us that Sie^^ (q) is transverse: 

q'Sa,{q)=0. (2.37) 

The array Siig{q) is therefore symmetric, transverse and traceless, and thus describes 2 physical 
tensor symmetries. Therefore, for n = we have in total 3 physical symmetries (1 dilation and 2 
anisotropic rescalings) . 



Here, we need only check that the nonhnear (field-independent) part of 6"/iio{q} is transverse. The field-dependent parts 
are guaranteed to be transverse by construction i.e. d'^Sj^Jg = (Section 2.2). 
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• 71 = 1: In this next simplest case, (2.32) requires 

MM = --Mui. (2.38) 

This gives 3 conditions on the 18 algebraically independent components of Afj^g^^. These include 
the 3 SCTs (2.8) of the scalar sector, with constant 6*: 

Mu% = hAto + K^ih - biSioir ■ (2.39) 

Under the SCTs, ( transforms non-linearly while ^a^ transforms linearly. In particular, (2.35) is 
automatically satisfied, for Mf^'^'^^ + Mf^^'^^ — ^SiegMf^"^ = 0. The non-linear shift in ( describes in 
real space the generation of a linear-gradient profile for the curvature perturbation. 

The remaining symmetries are "tensor" symmetries, 

Mill = Mm = (n = 1 tensor symmetries) , (2.40) 

under which (^ transforms linearly while 7 transforms non-linearly. The shift in 7 describes the 
generation of a linear-gradient tensor mode. The transversality condition (2.35) requires 

q' {Mi,,,, (q) + Me„ie, (q)) = . (2.41) 

This imposes 8 conditions on the 12 algebraically independent components satisfying (2.40), leaving 
us with 4 tensor symmetries. Therefore, for n = 1 we have in total 7 physical symmetries (3 SCTs 
and 4 tensor linear-gradient transformations). 

For n > 2, the counting of symmetries works as follows. The array Mn,...,^, being symmetric in 
its last n + 1 indices, starts out with | (n -|- 3) (n -|- 2) algebraically independent components. The trace 
constraint (2.32) gives ^{n -\- l)n conditions, while the transversality condition (2.35) imposes 3(2n -|- 1) 
relations on the coefficients. In total, we therefore have 

3 3 

-{n + 3){n + 2) {n + l)n - 3(2n -M) = 6 symmetries (2.42) 

at each order n > 2. These include 4 tensor symmetries, under which only 7 transforms non-linearly, plus 
2 "mixed" symmetries, under which both C and 7 transform non-linearly. 

2.5 Physical interpretation 

To shed light on the physical origin of these symmetries, consider expanding the metric about the origin 

00 ^ 
^^i = E^^*J^i-^"^''---^'"' (2.43) 

n=0 

where the array Hiji,...,^ is constant, and is symmetric both in its first 2 indices and its last n indices. At 
linear order, the metric is given by hij = (1 -|- 2Q)5ij + ■^ij. In particular, 7jj = hij — ^6ijh^^ is transverse, 
which implies the comoving gauge condition 

d^hij = ^dih', . (2.44) 
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This translates into a trace condition on the coefficients 

Him2-e.n = o-^««^2-£n • (2.45) 

(As a check, for the pure gauge configuration hij = di^j + dj^i = ^^^=0 ^(^ii4-^n +-^ii^i-€n)^^' ■■■x^", 
this trace condition reproduces (2.32).) 

As before, we restrict our attention to metric configurations which can be smoothly extended to a 
physical profile with suitable fall-off behavior at spatial infinity. In momentum space, the tensor profile is 



oo 



an 



^^M = E ^ {h,,,,...,^ - ^5.,/7«,,...,„) g^^^ ^^[ g^^ [i^Trf 6\q)j . (2.46) 

To enforce the finite-momentum transversality condition ct^ij{q) = 0, the Hij£-^...£^ must become q- 
dependent such that 

q' {Hi,,^...,M - \5,jHm,...iM)^ =0. (2.47) 

This is the generalization of (2.35). 

The counting works as follows. The array Hiji^...^^, being symmetric in its first 2 and last n 
indices, starts out with 3(n + l)(n-|-2) algebraically independent components. The trace condition (2.45) 
gives \{n + l)n conditions, while the transversality condition (2.47) imposes additional relations on the 
coefficients. At n = and n = 1, we find that the number of H coefficients matches the number of M's, 
respectively 3 and 7. This reflects the fact that hij and its first-derivatives can be made trivial at a point, 
which defines Riemann normal coordinates. At n = 2, i? has 12 components, while recall that M only 
has 6. This leaves us with 6 physical components, which matches the number of algebraically independent 
components of the Riemann tensor Rikjg. in 3 dimensions. ^'^ 

2.6 Noether charges 

The Noether charge associated with the symmetry transformation {^C^^Hj) is- 

g = ^y"d3:E(^{n^(x),5C(x)} + {n!/(x),57^,(^)}) , (2.48) 

where II^^ = 5C/5C, 11^ = 5C/5jij are the canonical momenta, and where the anti-commutator, denoted by 
{ , } , makes Q Hermitian in the quantum theory. The integration over an infinite volume may be divergent 
and the charge undefined, just as in theories with massless Goldstone bosons, but the commutator of Q 
with local fields is well-defined. In particular, it is straightforward to check with the above definition that 
the charge generates the symmetry transformation, [Q,C] = —iSC and [Q,'yij] = —i6^ij, as it should. In 



^'^For 71 > 3, we find fewer components than the most general coefficients in the Riemann normal expansion, which are 
given by derivatives of Riemann. For instance, a,t n — 3, H has 18 components, while M again has 6, leaving us with 12 
physical components. This is 3 fewer than the number of independent dmRikjt coefficients in the Riemann normal expansion. 
Although we have not checked this in detail, this must be because the metric profile we are considering is not the most general, 
but is restricted by the adiabatic transversality condition (2.47). This must translate in Riemann normal coordinates as 
adiabatic transversality constraints on derivatives of the Riemann tensor. 
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practice, we will regularize the charge by turning the integral over volume into a Fourier space quantity 
with the associated momentum understood to be taken to zero eventually. 

For the symmetry (2.27), the Noether charge Q takes the form 

- /■* df 
Q = Q + AQ —— , (2.49) 



H{t') 

where Q and AQ have no explicit time dependence. Q has explicit time dependence, but since it is 
a symmetry its total time derivative vanishes. We will be particularly interested in the part of these 
operators that generate the non-linear field transformations 

Qo = I d'^ 9,^,(f) Q5^^n^(x) + 2n!/(f)) = i I ^ q^,{-q) {^^5'm^{q) + 2Ii'^{q)^ ; 

AQo = Jd'x(^V^d,Cj{x))2Ul^^{x) = -iJ ^q^q,^,i-q)2Ui^{q). (2.50) 

For the diffeomorphisms ^^ in the Taylor expansion (2.31), the generators Qo ^^^ explicitly given by: 

qW = lim L^M,,„..,,„- -— -5^''^U^{q) + 2nf (g) , (2.51) 

g-^o n! dqij^---oqi„ \3 ^ J 

where we have used 11* ^ = 0. 

It is worth emphasizing that AQ is associated with the diffeomorphism ^* = V^^*, with f satisfying 
V^^j + did^^j/3 = (see (2.26) & (2.27)), implying the diffeomorphism ^* = V^^' is divergence free and 
is thus a tensor symmetry. As such, by the argument at the end of Sec. 2.3, the diffeomorphism receives 
no time-dependent correction and AQ by itself is a good conserved charge. 

3 Inflationary Consistency Relations as Ward Identities 

In this Section we derive the Ward identities associated with the non-linearly realized symmetries identified 
above. Analogously to the low-energy theorems for pions [38, 39], the Ward identities constrain the soft 
limit of various inflationary correlation functions. As particular cases of these, we will recover in Sec. 4 
the standard consistency relations [3-5], which determine the q^ and q behavior of the soft limits in terms 
of lower-order correlation functions. Beyond these, the Ward identities yield an infinite network of further 
consistency relations, which at each order partially constrain the g" behavior of correlation functions. 

The Ward identities are obtained by taking the in-in vacuum expectation value of the action of the 
charges 

{n\[Q,o]\n) = -i{n\so\n), (3.1) 

where |0) is the (Heisenberg picture) in- vacuum of the interacting theory, and 0{ki, . . . , fc^v) denotes an 
equal-time product of N scalar and tensor fields:^^ 

0{ku ...,kN) = OHku • . . , kM) ■ 07m+i,m+i,...,W.v(^m+1' ■■■^kN), (3.2) 

^^Since this is an equal-time product, the fields all commute with each other, hence their ordering is irrelevant. 
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where < M < iV, with O? = UaUCika,t) and O/^^,,-^^,,...,^,,^ = 0^=^+1 7^,,, (4, *)• Note that the 
tensor indices are arbitrary — a subset of these can be contracted among themselves or not. To avoid 
cluttering the notation, we will refrain from explicitly writing these indices unless necessary. We work 
within the Heisenberg picture unless stated otherwise. 

3.1 Time-independent identity 

Substituting the general expression (2.49) for Q, the Ward identity becomes 

- r* df' - - r* df 

{n\[Q,om + {n\[AQ,o]\n)J ^^^ = -i{n\60\n) - i{n\5AO\n) J -^^, (3.3) 

Although this explicitly depends on time, we now argue that the identity holds for the time- independent 
components Q and AQ separately. The argument is straightforward. First, consider a time- independent 
tensor symmetry ^T) ^-e., one under which 7,^ shifts non-linearly while C transforms linearly. By definition, 
this satisfies diS,^ = V'^^t ~ ^- ^^ discussed at the end of Sec. 2.3, such a diffeomorphism receives no 
time-dependent correction and represents an honest-to-goodness symmetry. The corresponding charge 
Q = Qt therefore satisfies 

{n\[QT,o]\n) = -i{n\6rrO\n). (3.4) 

In particular, since the generator AQ appearing in (3.3) is itself a tensor symmetry, as argued in Sec. 2.3, 
we have 

{n\[AQ,o]\n) = -i{n\6AO\n). (3.5) 

The time-dependent pieces in (3.3) therefore cancel out, leaving us with 

{n\[Q,o]\n) = -i{n\50\n). (3.6) 

Thus the Ward identity holds for any time-independent spatial diff^eomorphism ^^, subject to the gauge- 
preserving condition (2.7). In particular, it holds for each ^^''^' in the Taylor expansion (2.31). 

3.2 The left-hand side 

We begin with the left-hand side of (3.6): 

{n\[Q,o]\n) = -2iim{n\OQ\n) , (3.7) 

where we have used the fact that Q and O are hermitian. The operators are assumed to be evaluated at 
some (late) time t. For the most part, we suppress this t dependence, but its presence will be important 
in some of our arguments below. 

The in- vacuum \Q) of the full interacting theory (i.e., the Bunch-Davies vacuum) is related to the 
free vacuum |0) by: 

\n) = n{-oo)\o) , (3.8) 

where 

n{ti) = U^{ti,0)Uo{U,0), (3.9) 
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with U and Uq denoting respectively the fuU and free time evolution operators. This kind of statement 
should strictly speaking be understood in the context of a wave-packet (see [56] ) . In a similar way: 

n\-oo)Qn{-oo) = Qo, (3.10) 

where Qo is the free part of Q, and it generates transformations in (^ and 7 that are independent of 
the fields i.e. only the nonlinear transformations (2.50), since these are the symmetries of the action 
of quadratic fluctuations. This statement, that a Heisenberg operator sandwiched between 0(— 00) and 
its inverse equals its free part, is strictly true for a general operator only if it is evaluated in the far 
past [56]. However, if the operator happens to be independent of time, the time of evaluation becomes 
inconsequential. For our argument here, we treat Q as time independent, and apply (3.10) with Q 
evaluated at some time t. In reality this is not quite correct, and a full justification for our procedure is 
provided in Appendix B. 

Thus, we have 

Q\n) = n{-oo)Qo\0) . (3.11) 

The action of Qo on the free vacuum |0) can be worked out by using wave-functionals [57, 58]. Let us use 
ICo)7o) to denote eigenstates of the free Heisenberg operators Co and 70 (suppressing the indices on 70 to 
avoid clutter). Inserting a complete set of eigenstates, and using (2.50), we find:^^ 



QojO) = / i^Co^7oQo|Co,7o)(Co,7o|0) 



/ 



^^^^s1^^j{-Q} J DCoD^o (^Id^m^M + ^niiq)^ |Co,7o)(Co,7o|0) 
''' ,U-^1 1 DCoD,o\Co,,o) {-l''^Jch)-''h^) ^^-^°|0)- (^-12) 



(27r)3 
The free vacuum wavefunctional (Coi7o|0) takes the Gaussian form: 



(Co,7o|0) =AAexp 



-Co{k)D^{k)Co{-k) + j,,^{k)D^{k)^o''{-k) 



(3.13) 



(27r)3 \^2^"'^^'^^"' ' ' 4 
where J\f is an irrelevant normalization. The real part of the kernels D is fixed by the power spectra^^ 

Meanwhile, the imaginary part is related to the real part by the Schrodinger's equation on the wavefunc- 
tional [58], though we do not need its explicit form. Taking functional derivatives of the free vacuum 
wavefunctional we have 

QojO) = - I ^Q^U-Q) {^iDdlW'Coiq) + D,iqWiq?l |0) . (3.15) 



^^ Focusing on scalar perturbations for simplicity, the general relation in real space for the momentum acting on an 
arbitrary state \s) is Tl^(x)\s) = J D(^\(^)[—iS/5(^{x)]{(^\s). Fourier transforming into n^(g) — J d'^xll(^{x)e"''^ thus gives 
ndq) = jDC\(:}HS/Sa-mC\s}- similarly, {s\ndx) = / DC[i5/5C(f)](sK){CI implies {s\Udq) = / ^C[i<5/5C(-a)](s|C>(CI- 



^^See Appendix A for details on our conventions for the power spectra. 
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At this point we would like to repeat (in reverse) the same argument as in (3.11), so that we can 
pull the free Coio) and 70(g) through ri(— 00) to obtain the full Ciq) and 7(g). Recall that this argument 
strictly works only for (Heisenberg) operators in the far past, i.e., focusing on (^ for the moment: 

lim n^ti)C{q,ti)n{ti) = lim Co{q,ti) , (3.16) 

whereas what interests us are C and Co in the far future (which we denote by C{q) and Co(g) without any 
time argument). What saves us is that, ultimately, we are interested in the small-g limit, because the 
diffeomorphism in momentum space ^j (in 3.15) contains derivatives of the delta function d^{q). For any 
given tj, the following holds: 

mti) = C{q)[l + ^] , Co{q,U)=Co{q)[l + Ao] where lim A, Ao = , (3.17) 

i.e., C and Co approach a constant in the long- wavelength/super-horizon limit. Here, tj is the conformal 
time, to be distinguished from the proper time tj; the combination q\Ti\ serves to compare the wavelength 
with the horizon. ^^ The constancy of C, in the low momentum limit has been shown quantum- mechanically 
in [30, 31]. Plugging (3.17) into (3.16), we conclude that 

hm ri^{ti)C{q)n{U)=Co{q). (3.18) 

ti-5--00, (j|Ti|-5>0 

An analogous expression holds for the tensor mode 7 as well.^^ Thus, combining (3.11) and (3.15), we 
find 

^\^^ = -j (03^^^^(-^"^ i^^D^mKiq) + D,{q)^^\q)^ \^) , (3.19) 

where we have used ri(— oo)Co(q) = C(9)^(~oo), whose precise meaning should be understood as (3.18) 
(likewise for 7). One issue deserves more discussion: is (3.18) really adequate, given that we ultimately will 
be taking derivatives of C,{q) (hidden in ^j{—q)) before sending q to zero? Let us postpone the discussion 
to the end of this subsection. 

Substituting (3.19) into (3.7), and recognizing that i*j{q) = ^j{—Q)-, CHo) — C{~Q)^ iHq) = l{~'l): 
and D^ , D^ depend on only the magnitude of q, we arrive at 

{n\[Q,o]\n) = I ^q.^,i-q)(^^^{n\mo\n) + j^^ . (3.20) 

We have used the fact that O consists of a bunch of C's and 7's at the same time as Cio) and j{q), and 
thus they commute. For each term in the Taylor expansion (2.31), this gives 

(f7|[QW,0]|f]) = hm^^M,,„..., ^ (y^o(g-^O) + ,^U%0) . (3.21) 

q-^o n\ oqt^---oqi.^\P^{q) ^Pdl) J 

^^Note that t — )■ —00 corresponds to r — )> —00, while f — > 00 corresponds to r — >■ in de Sitter space. 
^^Note the somewhat intrincate limit implicit in (3.18): q should be sent to zero before ti is sent to —00. 
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Hence each Ward identity will constrain the g" behavior of correlation functions in the soft limit. This is 
the main result of this subsection. Let us close by returning to an issue raised earlier. From (3.21), it is 
evident that what we need is something stronger than (3.18) i.e., we need instead: 

lim Ot(t,)a," [f{q)Cm n{U) = 9," [f{q)Com , (3.22) 

where f{q) is some function of q, and dq^ represents some general n derivatives with respect to q (and 
likewise for 7). To justify this, we start from what we know to be true: 

hm n\u)dq''[f{q)aq,U)]niU)= lim dg"[f{q)Coiq,U)], (3.23) 

ti— >■— 00 ti— >— 00 

which follows from (3.16). Substituting (3.17) into dq'^[f{q)Ciq,ti)], we have 

d.^'immU)] ~ dq^[f{q)C{q}] + dq^[f{q)m]^ + dq--'[f{q)Cmdg^ 

+dq--V{q)m]d,'^ + ..., (3.24) 

where we have been cavalier about numerical coefficients. The important point is that dq'^~'^[f {q)(^{q)] ~ 
q~"^dq"'[f{q)({q)], and thus as long as 

lim g™a„™A = for m > (3.25) 

(and likewise for Aq), all the corrections from A, Aq and their derivatives drop out in the long- wavelength 
limit, and (3.22) is established. Therefore, the precise long-wavelength-constancy requirement on (" is 
(3.17) supplemented by (3.25). It is worth noting that (3.25) is obeyed by many different possible A's. 
For instance, 1 + A = (l + zgrj)e~*'^'^* (i.e., Hankel function) gives A ~ (grj)^/2 + ... in the long-wavelength 
limit, while 1 -|- A = e*''^' gives A ~ iqTi + .... Both satisfy (3.25). The former describes the behavior 
in standard inflationary models, while the latter arise in more unusual scenarios including ones with no 
expansion [41]. It is important to emphasize that whether A goes like q or q^ in the soft limit has no 
bearing on which of the consistency relations {i.e., which n in (3.21)) is satisfied. As long as (3.25) holds, 
the consistency relations are expected to hold for all n. 

3.3 The right-hand side 

At each order in n, the variation of O can be split into a part that includes the non- linear shifts in the 
fields and a part that includes the linear transformations: 

5-WO = 5ll_,^0 + 4">. (3.26) 

As shown in Appendix C, the non- linear part contributes to disconnected diagrams in the Ward identities 
and therefore drops out when considering connected correlation functions. For instance, from (2.33) the 
non-linear variation of C, schematically gives contributions of the form 






(Ci-lin.C^> ~ Ijm 1^^\?+ kaHCikl) ■ ■ ■ Cika-l)aka+l) ■ ■ ■ C(fcM)0^(/^M+l, • • • , %)) 
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describing a 2-point function (or more precisely, {Cio) C (ka)) / P(iQ)) disconnected from an A^ — 1 point 
function, which is precisely canceled by a disconnected term from the left-hand side of the Ward identities. 
Focusing on connected correlators, substitution of the field transformations (2.33) therefore gives 

^ an 

_ J2 r^°^°^°(fc,) ^^^ _ _ .^^^ (0^(fcl, . . . , 4-1,4+1, • . . kMH.jAka)0^kM+l, ..., kN))c 

^ an ) 

- E r-\>^4) , ,^ , (0^(4, . . . , 4/)o7^,,,,,,,,...,,,,„...,,,,(4.+i, • • • , a?^))c + 

6=M + 1 ^1 ^n J 

(3.27) 
Combining with (3.21), we obtain the connected version of the Ward identities: 

HmM,,„..., ^\ (^AY''{q)Oiku • • • , %))c + ,^AmO{ku • • • , %))c ) 

a™ ki d 



iV / r^n I i Cin+l 



_ ^ r^o..>(^j______^0^(4, . . . , 4-1, 4+1, . . . 4/)7wa4)o^(^M+i, . . . , ^^)>c 

a=l "l *^n 

- E r\>^n4) , . , (0^(4, . . . , ^^m)oL,,,.,,....,.,4,......(^Vi, . . . , A?^))c + . . . 

(3.28) 

The ellipses indicate terms on the right-hand side that are higher-order in the fields. One further simpli- 
fication is possible. As shown in Appendix D, the momentum-conserving delta functions implicit on both 
sides of (3.28) can be canceled, leaving us with an identity for "primed" correlation functions, defined by 
removing the delta function [47] 

{0{q, ki,...,kN)) = {27Tfd''{P){0{q, 4, . . . , kN))' , (3.29) 

where P = q + ki + . . . kjsf. The primed correlator is thus defined on shell. It is a function of only A'^ 
momenta, which can be chosen to be fci, . . . , k^. As shown in Appendix D, the Ward identities in terms 
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of primed correlators are given by 

(9" / 1 ^ ^ n ' ^ - 

_ ^ y^o...-(fcJ ^^^__ .^^^ (0^(fci, . . . , 4-1, 4+1, . . . kMH^,AkaW{kM+l, . . . , fciv)); 
a=l ^1 "n, 

N ^n ) 

b=A/+l ^1 ^n j 

+ . . . (3.30) 

Apart from the 6no term on the second hne, this takes the same form as the unprimed identity (3.28). 
The general Ward identities (3.30) relate N + 1-point correlation functions with a g = scalar or tensor 
insertion to the symmetry transformations of A^-point correlation functions. 

3.4 Component Ward identities 

The Ward identities (3.30) involve the coefficient arrays Mj^(,,..£^(g). In order for these to arise as the g — )• 
limit of correlation functions involving physical modes, the Ma^^^j^ parameters of the transformations 
must acquire g-dependence, as discussed in Sec. 2.4. Specifically, the regulated array Mj£g..,^^(g) must 
satisfy the transversality condition (2.35). For each choice of M's satisfying these conditions as well as 
the symmetry property (2.32), we can plug into (3.30) and obtain a Ward identity. 

Alternatively, we can remove the Mj^g,..£,^(g) coefficients from (3.30) by projecting the indices con- 
tracted with the M's onto an appropriate subspace. We must not only ensure that the resulting identities 
are symmetric in (£o, • • ■ , ^n) and consistent with the trace condition (2.32), but that they are also trans- 
verse in the sense of (2.35). This can be achieved by introducing operators Pie.Q...e.^jmo...m„{(l) with the 
following properties: 

1- Piio...i„jmo...m„ is symmetric in the {Iq, • • • , ^n) indices and in the (mo . . . rUn) indices. 

2. Piio...injmo...m„ is symmetric under the interchange of sets of indices: Piio...i„jmo...m„ = Pjmo... mnilo ■■■£«■ 

3. For n > 1, PiiQ...i„jmo...m„ obeys the trace condition (2.32): 

Pim2...i„jmo...m„ = -^Pm£2-enjmo...m„ (n > I) . (3.31) 

4. Piio...i„jmo...m„ satisfies the transverse condition (2.35): 

9* {Piioii...i,ijmQ...mniQ) + Peoiii-i«jmo...mn{Q) " o '^««0-f'«^l-Mmo...m„ (^) j = ^ ' (3.32) 

In other words, P has the same properties as M under either sets of indices and is symmetric under 
the interchange of sets of indices. In Appendix E, we will explain how to systematically construct these 
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physical operators and give explicit expressions for the first few values of n. (As noted in the Appendix, 
the various projectors obtained using this method will in general form an over-complete set. This means 
the identities below will not all be independent, though none will be missing.) 

In terms of these operators, the Ward identities (3.30) become 






dk'^,---dk'^^ N n + ldk^^---dk^J 



^^ an 

_^T^--oW.(^j^ {0<^(k^,...,K-,,K+l,.■■kMH^aka)0^kM+l,...,kN)yc 



(3.33) 
These component Ward identities are the main result of this paper. A few remarks are in order: 

• At each n, the identities (3.33) constrains the g" behavior of the soft limits. At lowest order, n = 0, 1, 
the q^ and q behavior is completely fixed [3, 43]. At higher order, n > 2, the Ward identities only 
constrain part of the A^ + 1 correlator in the soft limit. 

• The primed correlation functions are defined in (3.29) as on-shell correlators. Correlators on the 
left-hand side, such as {({q)0{ki, . . . , k]^))' , are functions of N momenta. Without loss of generality, 
these can be chosen to be q,ki,..., kjy^i, in which case one should make the replacement kjy = 
—q — ki — . . . — /cjv_i before differentiating with respect to q. Similarly, the primed correlator 
{0{ki, . . . , k]\})y on the right-hand side is a function of A^ — 1 momenta. Choosing these to be 
ki, . . . , kfyf-i, one should make the replacement k]sf = —ki — . . . — k^-i before differentiating with 
respect to ka, where a = 1, . . . , A^ — 1. 

• The ellipses denote corrections that are higher-order in the fields, as discussed below (2.33). (As ar- 
gued in Sec. 2.2, however, the dilation symmetry is exceptional, and the field transformations (2.16) 
are exact in this case.) 

• The Ward identities (3.33) hold independently of whether the hard modes with momenta ki, . . . , k^ 
are inside or outside the horizon. In particular, following [37] it would be interesting to generalize 
the analysis to include correlation functions with time or spatial derivatives acting on the short 
modes, as this can be useful for loop calculations [59-63]. 

• Note that even if the operator O contains no 7, the right hand side does involve 7 through the third 
line of (3.33), which replaces successively each C, by 7. However, the fourth line would vanish in 
such a case. 

• The Ward identities also have implications for correlation functions with a soft internal line, where 
the sum q = —{ki + . . . /cjv/) of M of external momenta approaches zero. In this limit, the correlation 
function is dominated by the exchange of a soft scalar or tensor mode, and the answer factorizes 
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into the product of A^ — M + 1- and M + 1-point correlators in the soft hniit [64, 65]. For instance, 
for a product of scalar modes [43] , 

Um(0^(fci,...,fc;v)>'e = lim ( (C(glO^(^i, • • • 4/))'cTrTT(C(-g10^(4m, • • • M>c 

+ ^(7«(g-)0^(^i,...^,,));-^(7«(-g-)0'^(4,+i,...fc^))M .(3.34) 

The identities (3.33) can be applied to the soft limits of the individual correlators in the products, 
and thus constrain the form of correlation functions with soft internal lines. 

In the following Sections, we will study special cases of these identities. Specifically, we will show in Sec. 4 
that the n = identities reproduce Maldacena's original consistency relations for scalars and tensors [3]. 
We will also show that the n = 1 identities reproduce the conformal consistency relation and the linear- 
gradient tensor consistency relation derived recently in [43]. In Sec. 5, we will show in detail that the 
n = 2 identities lead to novel consistency relations, and check their validity with the 3-point graviton 
correlation function in slow-roll inflation. 



4 Recovering Known Consistency Relations 

In this Section we will recover, as particular cases of our general Ward identities (3.33), the known 
consistency relations for single-field inflation. Specifically, for n = we will recover Maldacena's original 
dilation consistency relation and anisotropic scaling relation [3]. For n = 1, we will reproduce the linear- 
gradient consistency relations derived recently in [43]. 

4.1 Dilation consistency relation 

For a dilation, the transverse condition (2.35) is trivially satisfied since tensors transform linearly. We 
can simply go back to (3.30) and substitute Mj^q = XSn^ to obtain^^ 



hm-l^(C(g10(^i,...,M)'e 



( 3(iV - 1) + f; 4 . J^ ) {0{h, ..., kM))', , (4.1) 

V a=l dka) 



which is the well-known consistency relation of inflation [3~5] . Note that this identity receives no higher- 
order corrections in the fields, since the dilation field transformations (2.16) are exact. 

4.2 Anisotropic scaling consistency relation 

For the anisotropic scaling, we substitute in (3.33) Mj^^ = 5j^q with 5*^ = and q'^Sa^^q) = 0. For 
concreteness, let us specialize to the case where the hard modes are scalars: = 0^ = Y[a=i C(^a)- In 



^® Alternatively, one can work directly with (3.33) by substituting Pugjmo ~ h^Uo^i^^oj which has all the desired properties. 
Tracing the resulting identity over {i,£o) yields (4.1). 
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this case, the Ward identity reduces to 

1 ^ i 8 

- -{0^{kl, . . . ,ka^l,ka+U . . . ,kN)ljmo{ka))'c \ 

+ ... (4.2) 

where we have used '^^^^^°'-''^{ka)^iaja{ka) = \ljmo{ka). The required operator Pji j^ has the same 
permutation properties as M for each pair of indices, that is, it is symmetric and traceless in (i,^o) and 
{j,mQ). Moreover, it is symmetric under the exchange of pairs of indices {i,io) -H- (j, w-q), as well as 
transverse q^P^ ,„j = 0. The operator with these properties and suitably normalized can be readily 
inferred: 

Pieojmo ~ PijPlomo + PimoPjlo — PiioPjmo , (4.3j 

where Pj£q = Sug — qiqe^ is the transverse projector. (The operator P^ ■ appears in the completeness 
relation for the polarization tensors — see (A. 7).) Since 7jmo is itself transverse and traceless, we have 
P^iojmoiQh'""'iQ) = 27i£o(€), and (4.2) simplifies to 

1 -- 1 ^fa-- 

- -(0^(A?i,...,4-i,4+i,...,A?iv)7jmo(4)>c> + ... (4.4) 

Since this relation is usually expressed in the helicity basis, reviewed in Appendix A, we can project both 
sides with the polarization tensor ef^ (?). Using the orthonormality condition (A. 6), we obtain 

- -{0'^{ki,...,ka-i,ka+i,...,kNhieuika))'c\ + ■■■ (4.5) 

For N = 2, the last line vanishes (since (C7) = 0), and the result agrees with [3]. For general A^, the 
above agrees with [43] to lowest order in the tensors, that is, as long as the last line is neglected. 

4.3 Linear-gradient consistency relations 

It has been argued recently that the order q behavior of soft correlators is also fixed by consistency 
relations [43]. We will show that these relations follow from the n = 1 Ward identities. To make contact 
with [43], we again specialize to a product of scalars for the hard modes: O = O^ . 
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Recall that at n = 1 we have a total of 7 symmetries, consisting of 3 SCTs and 4 tensor symmetries. 
Starting with the SCTs, we substitute MfpJ = bi-^Sng + be^Sa-^ — biSi^i^ in (3.30) to obtain 

q^o dq" yPc_{q) J ^^^\dK dkidki dkidkl) 

+ .... (4.6) 

(Note that the T terms in (3.30) vanish identically in this case.) This agrees with the conformal consistency 
relation — see Eq. (54) of [43]. Although it was originally believed that the SCTs were restricted to the 
scalar sector only [43, 44], we have now generalized these transformations to include tensors. In particular, 
the tensor corrections encoded in the ellipses in (4.6) can be computed explicitly using the expressions 
given in Sec. 2.2. 

For the tensor symmetries, described by fully traceless Mj^^^^, the Ward identity (3.33) gives 



= -Pli,,^,mM) E \ \^^ dk^\k^ ^^'^^^' • • • ' ^^^^'^ 

_r^mo^.j^^ka)^{0<{k^, ..., 4-1, K+l, ..., kNH,,{ka)t [+••■, (4.7) 

where PjiQi^jmami ^^ symmetric in (^0,^1) and (mo, mi), fully traceless in {i,(.Q,(.i) and (j, mo,mi), sym- 
metric under the interchange {i,io,ii) -H- (j, mo,m,i), and transverse: 

^ ( PitQlijmormW + -ffoi^ljmomi (9) ) ~ ^ ' ^^'"^ 

As described in Appendix E, we may construct three linearly independent operators satisfying these 
conditions: 

jT (1) / ~N _ pT pT 



p^ (^) (fi\ ^ pi pi 

pT (2) ,~x ^ pT pT , pT pT -^rpT pT , pT pT ^ 

ilo£ijmom-i\^J ikjmi-^ lolimok '^ ^ ikmoj^ £oiimik o V iokmomi-^ Hijk '^ likmomi-^ f-o^jk) 

pT (3) /^N _ ^j t -~l^ pT I -£0 pT _ -J pT \ I -mo / -£i pT , -£o pT _ -i pT 

where P'^i^j^n i^ defined in (4.3). To cast (4.7) in the helicity basis, we focus on P^^ ^ ■ ^ and con- 
tract both sides with qi^^ef^ (q). Using the orthonormality (A. 6) and completeness (A. 7) relations of the 
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polarization tensors, the result is 

X (0^(fci, . . . , 4-1, 4+1, . • • , kNHeoiQYc > + ■■■ (4.10) 
This is consistent with [43] to leading order in the tensors. 

5 Example of Novel Consistency Relation 

The first novel consistency relation arises for n = 2, which constrains the g^ behavior of correlation 
functions in the soft limit. Unlike the n = and n = 1 consistency relations of Sec. 4, which fully 
constrain the q^ and q^ behavior of the correlators, the n = 2 Ward identities only partially constrain the 
q^ behavior of the correlators in the soft limit. 

To illustrate this, let us focus on the n = 2 "tensor" symmetries, for which Mn^i-^i^ is fully traceless. 
The corresponding P'j£f^£-^£,jjmnmim2 ^^ symmetric in the (^0,^1,^2) indices and in the {1120,1711,1712) indices, 
traceless in the (i, ^0,^1, ^2) indices and in the {j,mQ,mi,m2) indices, and symmetric under the exchange 
(i, £0,^1,^2) ■^ (i, "io,"ii, "i2), and transverse in the sense of (3.32). See Appendix E for a description 
of how to arrive at operators with these properties. Furthermore, to allow comparison with existing 
computations, we specialize to the case of A^ = 2 hard scalar modes, 0{ki, k2) = Ct C^ . In this case, the 
T terms are absent since (C7) = 0. 

With these assumptions, the Ward identity becomes 

1- pT / -N, d I I j„jQ,^ r \' 1 — — P""" ( ") \^ ^^ki^k2' 

(5.1) 

Note that we have neglected the "..." corrections which are higher-order in the fields. In slow-roll inflation, 
these corrections are higher-order in the slow-roll parameters and hence can be consistently neglected 
to leading-order in slow-roll. We will check this identity using the three-point function computed by 
Maldacena [3]^'^ 

j4j(...(.ic,c,)-<.-«)^^i*r(-.-.^^^±^t±^.^) . (-, 

where K = q + ki + k2, and -fiXiimn(^) ^^ defined in (4.3). Meanwhile, the scalar two-point function is 



Aekf 



iCkSkJ' = T7U- (5-3) 



^^To translate from Maldacena's expression in the helicity basis, we multiplied both sides by tjrno (9), summed over helicities 
and used the completeness relation (A. 7). 
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Both (5.2) and (5.3) are valid to leading order in slow-roll parameters. Before verifying (5.1), we quickly 
check that the lowest-order consistency relations are satisfied. 

• Anisotropic scaling consistency relation (n = 0): In the limit g — )■ 0, we have /c2 — s- —ki, and the 
left-hand side of (4.4) gives 

S P^^i-r'^mshY = f5.,™.(«)^5*i^r . (5.4) 

Differentiating (5.3), the right-hand side of (4.4) becomes 

which agrees with (5.4). Thus the n = Ward identity (4.4) is satisfied. 

• Linear-gradient consistency relation (n = 1): To compute the g-derivative of the three-point func- 
tion, we must let k2 = —ki — g and work consistently to linear order in q: 

+ (higher — order in q) . (5.6) 

Next we contract this with P^ e jm m • ^^ doing so, we will use the following identities, which 
straightforwardly follow from the properties of P^ ^ ■ ^ and the explicit form (4.3) for Pj^ j^f 

pT jmomi d pT _ r, . pT jmo pT _pT i pT /Tr 7\ 

Applying this operator, the left-hand side of (4.7) becomes 

imj f>i,,,„„„, iA) ^ (^ (7-- mi, ■■■ cj„ »::) = - y ■ i^ ■ *i^r^-r" ^i„...„.„. (« . 

(5,8) 
Meanwhile, on the right-hand side of the identity (4.7), we have 



1 ^ • 3^ 15 iJ^ - • ' - 

2^if.QhjmomM)Z_^^iQ^a QJ^a ^^kSk2^' ^ "Y ' 4g^ " ^l^T'K'^ PiloiijmoraM) ' (^-9) 



where we have used the fact that -PJn^^jmnmi ^^ traceless. This agrees with (5.8), which verifies the 
n = 1 tensor Ward identity (4.7). 

Novel tensor consistency relation (n = 2): We now verify that Maldecena's 3-point function (5.2) 
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satisfies the n = 2 Ward identity (5.1). On the left-hand side, we find 



1712 



2ki \ 5g™i 5g 

+ 4 A:i 1^ 1 ag-i ^1 ag™2 






^ ;^2 -Pjmofc^(^)0mim2 + „, 2 1 1 1 1 ^jmokliQ) 



+ (higher — order in g) . (5.10) 



When contracting this with Pj^f^i-^i2Jmomim2^ ^^ ^^^^ ^^^ ^^^ identities 



if.0e.-it2 inim,2 Q„inxQ„m2 ' i^o^ih mim2 Qgrrii 



pT jmo pT _ pT , pT /c;ii\ 

which follow from the properties of Pj e e jm m m ^^'^ ^^^ explicit expression (4.3) for Pj^ ^^. It 
turns out that only the last term in (5.10) gives a non- vanishing contribution: 

52 / I \ i72 35 _ ^ ^ ^ 

(5.12) 

Meanwhile, the right-hand side of the Ward identity (5.1) becomes 






/ ^ Q /3l.a /5i[,a flLa Aci-S /,2 

a=l '-"^mo'-"^mi'~"^m2 ^^'''1 '^1 

(5.13) 
Contracting with P^^ e e km m m. ^ ^^'^ using the fact that this operator is traceless, we obtain 



P 



T 
ieo£i£2Jmom 



^a dHC^^^C^J, _ H' 35^, 



a=l 1 ^ 



This matches (5.12), which verifies the n = 2 tensor Ward identity. 

Although we have focused for concreteness on the tensor symmetries in this check, there are of 
course 2 additional "mixed" symmetries at n = 2, under which both scalar and tensor modes shift non- 
linearly. It will be very interesting to study this identity in detail, particularly in models with reduced 
sound speed. 
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6 Conclusion 

In this paper we derived an infinite number of consistency relations constraining at each n > the g" 
behavior of cosmological correlation functions in the soft limit g — )• 0. We showed how they arise as the 
Ward identities for an infinite set of global symmetries, which are non- linearly realized on the perturba- 
tions. As the lowest-order identities (n = 0, 1), we recovered Maldacena's original consistency relations for 
scalars and tensors [3], as well as the recently-discovered linear-gradient consistency relations [43]. The 
higher-order (n > 2) identities are new, and we checked as a particular example that the n = 2 "tensor" 
identity is satisfied by known correlation functions in slow-roll inflation. Our general Ward identities 
hold whether the hard modes are inside or outside the horizon, and also have implications for correlation 
functions with soft internal lines. 

There are many directions that would be interesting to pursue: 

• To offer further insights on the higher-order consistency relations, it would be interesting to red- 
erive these new identities using standard background-wave arguments. This may be technically 
challenging already for the n = 2 identities, but would give an independent check of the Ward 
identities. 

• It would be enlightening to apply the methods developed here to study correlation functions with 
multiple soft external lines. In the case of pions, it is well-known that the double-soft limit probes 
the non-abelian nature of the broken algebra [39]. (See Sec. 4.1 of [66] for a recent discussion of 
double-soft pion theorems.) In the present context, focusing on the scalar sector with symmetry 
breaking pattern (1.3), correlation functions with two C's taken to be soft should similarly instruct 
us about the underlying broken de Sitter isometries. To zeroth-order in the soft momenta, this has 
recently been discussed in [37]. However, the non-abelian nature of the conformal algebra should 
show up at linear order in the soft momenta. This is currently in progress [67]. 

• Our derivation of the consistency relations from Ward identities crucially relied on the standard 
Bunch-Davies vacuum. Specifically, this assumption was made in Sec. 3.2 in choosing the vacuum 
wavefunctional. However, the derivation straightforwardly generalizes to arbitrary initial conditions, 
by substituting a modified wavefunctional, such as those of interest recently [68-73], in lieu of (3.13). 
The resulting identities would yield modified consistency relations with which to test non-standard 
initial conditions [74, 75]. 

Our derivation of consistency relations from Ward identities is very general and can be readily 
applied to other symmetry breaking patterns, including non-infiationary ones. The conformal mecha- 
nism [76-82], in particular, relies the spontaneous breaking of conformal invariance on approximate flat 
space, with symmetry breaking pattern 50(4, 2) — )• 50(4,1). The non-linear realization of 50(4,2) im- 
plies novel consistency relations, which have been derived recently using the background- wave method [83]. 
It should be straightforward to reproduce these relations from Ward identities. 
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A Conventions 

In this Appendix we briefly summarize our conventions. 

A.l Power Spectra 

The curvature perturbation is expanded in terms of Fourier modes as 

(^C(fe,t)e-^'=-^ (A.l) 

The scalar power spectrum P^(A;,t) is defined as 



{n\C{k)C{k')\n) = {2TTy6'\k + k')p^{k, t) . (a.2) 

We also define the free-theory power spectrum, P9{k,t), in terms of the 2-point function of interaction- 
picture fields with respect to the free-theory vacuum state: 

(0|Co(^)Co(^')|0) = {27Tf6\k + k')P^{k,t). (A.3) 

In slow-roll inflation, for instance, P9 = H'^/AMp^ek^. 
Similarly, we expand tensor modes as 

/ri^k 
^^7.,(fc,i)e-*^ (A.4) 

In the helicity basis, we follow [3] and write 

j,,{k,t) = Y,4,ikh'{k,t), (A.5) 

s=± 

where the polarization tensors e^- are transverse and traceless, ^*el- = efj = 0. These satisfy the orthonor- 
mality condition, 

e!jik)et-Ck) = 25ss', (A.6) 

and the completeness relation, 

J2^!j{k)eU{k) = P^'.keCk) , (A.7) 

s 
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where the projector P^i^^ is defined in (4.3). 
The tensor power spectrum is defined as 

{n\-f'{k)-f''{k')\n) = {2TTf6^{k + k')P^ik,t)5ss', (A.8) 

or, equivalently in terms of "jij, 

{n\-fij{k)j,j{k')\n) = {2^f5\k + k')AP^{k,t) . (A.9) 

The free-theory tensor power spectrum is similarly given by 

(Vt \loij{k)lOij{k')\ n) = {27rf6^{k + k')AP^{k,t) . (A.IO) 

In slow-roll inflation, for instance, P^ = H'^/Mp^k^. 

A. 2 In-in formalism 

We briefly review the in-in formalism. Our goal is to calculate the correlator {i}\0{t)\Q) for some operator 
O, where so far all quantities are in the Heisenberg picture. The operator O satisfies the Heisenberg 
equations of motion, O = i[H{t)0{t)], where H{t) is the full time-dependent Hamiltonian. The state \il.) 
is the state of the system, and does not depend on time in the Heisenberg picture. (This means that it 
is not an eigenstate of H{t) at all times, since these eigenstates would in general depend on time). 

As usual we split the Schrodinger-picture Hamiltonian into a free part and an interaction part, 
H{t) = HQ{t) + V{t), and move to the interaction picture with respect to this split. We have 

Mt))i = U^o{t,0)U{t,0)\n); 
Oiit) = U^Q{t,0)U{t,0)0{t)Uo{t,0)U\t,0). (A.ll) 

where U and Uq are the time evolution operators for the full and free Hamiltonian, respectively. Interaction- 
picture operators evolve according to the free Hamiltonian: Oj = i[{HQ)i{t),Oi{t)]. Interaction-picture 
states evolve according to a Schrodinger equation with Vj{t) playing the role of a Hamiltonian, 

i^^\nit))i = Vi{t)\nit))j, (A.12) 

with general solution 

\n{t))j = Ui{t,u)\n{ti))j; 

, (Texp\-iJldt'Vi{t')} t>ti, 

I Texp I i/;Mt' !//(*')} t<ti. 

The expectation value is the same in any picture: {i^\0{t)\0,) = j{il{t)\0[{t)\Q{t))i. Using the 
above expressions for the time evolution, we can write our correlator in terms of the state at the initial 
time ti, 

(fi|C(t)|Q) =/(f^(t)|0/(t)|fi(t))/ = /(0(t,)|C/|(t,ti)0,(t)C/,(t,ti)|f^(ti))7. (A.14) 
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The statement that \Q) is the in-state is the statement that in the far past it was the Fock vacuum |0) of 
the free field creation/annihilation operators, that is, 

hm \n{ti))i = \0). (A.15) 

To see that this coincides with the usual definition of the in-state, we can write the same expression as 

|fi) = lim U\ti,0)Uo{ti,0)\0) = n{-oo)\0), (A.16) 

which matches the usual definition. With this, we obtain the in-in form of the correlator (A. 14): 

{n\0{t)\n)= lim {0\UJ{t,ti)Oi{t)Uiit,U)\0) . (A.17) 

B Charges and Power Spectrum 

Here we justify a number of steps made in Sec. 3.2. Recall from Sec. 2.6 that the true conserved charges 
are AQ, and Q = Q + f{t)AQ where f{t) = f dt' /H{t'). As such, the arguments made in Sec. 3.2, 
especially (3.11), can be rigorously applied to them to obtain the analogs of (3.19): 

+f{t) I ^m'Cj{-q)D,{q)f^{q)\n) . (B.l) 

On the other hand, Q\Q,) = Q\Q) + f{t)AQ\Q). Thus, we are led to the conclusion that 

Q\^) = -j ^,^U-^ {^^D^mKiq) + D,{q)Y'{q?j l^) , (B.2) 

justifying (3.19). No assumption has been made about Q's time-dependence (or lack thereof) . 

Another subtlety glossed over in Sec. 3.2 is that we did not distinguish the C or 7 power spectrum 
of the full theory from that of the free theory. The reason has to do with the soft limit. In the g — )• 
limit, C{q) and Jijiq) are time- independent. This once again allows us to apply the analog of (3.11) at 
any time t we wish: 

cmn) = n{-oo)Co{q)\o) , Ji,{q)\n) = n{-oc)Jo^M\o)■ (B-3) 

We thus have 

{n\C{-q)Cmn) = {0\Co{-q)Co{q)\0) , {nhij{-q)^,j{q)\n) = {Ohoiji-q)lOijmO) (B.4) 

in the soft limit. 
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C Connected Ward Identities 

The connected correlators (/ci/c2 • • • k^Y ^"^^ defined as 

{k,k2 ■ ■ ■ k^r = E (1^1 - 1)' (-1)'"'"' n ( n ^'^ ) ' (c.i) 



BeiT \aeB 



where vr is a partition of 1, • • • ,N, \it\ is the number of blocks in the partition, and B represents a block 
of the partition. Here we abuse the notation a bit, and use fei, ..., A; at as stand-in for scalar or tensor fields 
at the corresponding momenta. The full correlators (^1^2 • • • k]\f) can then be written in terms of the 
connected correlators {kik2- • • kjy)'^ as 

{kik2---kN) = Y.ll/Ylka\ . (C.2) 

TT Ben \aeB I 

In our case we have (fc) = {kY = 0. Thus we can write 

{kik2---kN)=Y,Wl\{ka\ . (C.3) 

7r' Be-ir' \aeB I 

where the sum is only over partitions vr' where each block has at least two elements. The first few instances 



are 



{kik2) = {kik2Y, (C.4) 

{kik2ks) = {kik2ksY, 

{kik2k3ki) = {kik2kskiY + {kik2Y{kzkAY + {kiksY{k2k4Y + {kikiY{k2k3Y , 

{kik2k3k/ikrYi = (A;iA;2fc3^4fe5)'^ + [(^i^2)'^(^3^4^5)'^ + 9 more permutations] , 

(^1/02^3 ^4 ^5 fee) = {kik2k3kik^kQY + [{kik2Y {k^k^^k^k^Y + 14 more permutations] 
+ [{kik2k3Y {kik^k^Y + 19 more permutations] , 



Each connected correlator comes with an overall delta function, 

{kik2 ■ ■ ■ k^Y = {2TTYS\ki + fc2 + • • • + kN){kik2 • • • %) '. (C.5) 

The {kik2 ■ ■ ■ k^)' are defined on shell. The power spectrum is P{k) = {k, —k)'. 

Our general Ward identity relation for A^ > 1 is in terms of the full correlators, 

/I \ ^ 

N 

-Y,Da{ki---kN), (C.6) 



a=l 
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where Da and Da are differential operators that depend on and act only on the a**^ momentum. (For 
example, for the dilations we have Da = 3 + fca • g^-) For A'" = 1, the second line of the RHS does not 
contribute, the first line of the RHS gives (27r)^Pi(5^(/i:i) and the relation is satisfied. 

What we would like to prove is the following relation for connected correlators, for N > 2. 

N 

(C.7) 



lim Vg (-^{q fci • • • kN)") = - V Da{ki ■ ■ ■ k^Y . 
g^O \P[q) J ^ 



For N = 2, the first line of the RHS of the relation (C.6) vanishes, and there is no distinction in the 
remaining correlators between connected and full, and so (C.7) is true for N = 2. This is the first step in 
a proof by induction that (C.7) is true for N >2. 

For the induction step, assume that (C.7) is true for A'^ — 1. We start by breaking up the correlator 
{qki- ■ ■ kj^) on the left-hand side of (C.6) as follows 

N 
{qki-'-kiy) = {qki--- k^Y + ^^(^ ^a) {ki--- ka-ika+i • • • A^at) . 

a=l 
*' \ aeB I Be{TT'\B) \ae-B / 

Some explanation is in order: we have broken {qki- ■ ■ k^) into its connected pieces. The first term is 
the fully connected piece {qki- ■ ■ k^Y- The sum in the first line of the right-hand side is over all the 
partitions in which the soft momentum q appears in a two point factor. The rest of the partitions are in 
the second line of the right-hand side: here vf' is the set of all non-trivial^^ marked partitions of 1, 2, • • • , iV 
where each block has at least two elements. A marked partition is a partition in which one of the blocks 
of the partition, the one denoted by B, is distinguished. In this case, the special block is the one whose 
connected factor contains the soft momentum q. We denote by Tt'\B the set of blocks of the marked 
partition which are not marked. For example, with A'^ = 3 there are no terms in the sum over marked 
partitions (since in this case there are no partitions where all the blocks have at least two elements), and 
for A^ = 4, the sum reads 

{kik2){qk^kiY + {kikz){qk2kiY + {kiki){qk2kzY 

+ {k3k^){qkik2Y + {k2k4){qkik3Y + {k2k3){qkik^Y ■ (C.9) 

Looking at the summation term in the first line of the right-hand side of (C.8), using {q ka) = {2'KY6^{q + 
ka)P{q), and inserting into the left-hand side of (C.6), we precisely reproduce the first line of the right- 
hand side of (C.6). 

Thus, to complete the induction step, we need only show, under the assumption that (C.7) is true 
for < A^ — 1, that the second line of the right-hand side of (C.8), inserted into the left-hand side of (C.6), 
precisely reproduces the second line of the right-hand side of (C.6) (minus the connected part). The 



^*A partition is non-trivial if it does not contain the trivial partition consisting of all the elements, which is captured 
already by the connected part (gfci • • • fcjv)'^. 
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second line of the right-hand side of (C.8), inserted into the left-hand side of (C.6) and using (C.7) yields 



7f' \ ^^^ \ ae5 / / B£iTT'/B) \aeB 



E 






n n^-. 

B&{tt'/B) \aeB 



(CIO) 
Next we break up the correlator in the second line of the right-hand side of (C.6) into its connected pieces, 
where the sum is only over non-trivial partitions vr' where each block has at least two elements, 



{k^k2 ■■■kN) = {klk2 ■ ■ ■ k^Y + Y.X{\X{^- 



(C.ll) 



The operator Ylia=i ^o. can be broken up over a given partition vr' as YliBeir' X^aeB -^«' ^"^^ when applied 
to Yin' riBe-TT'dlae-B ^aY-: ^^ generate the sum over marked partitions. 



E 



E"«(n*^" 

aeB \aeB I 



n n*. . 

Be(n'/B) \a&B I 



(C.12) 



where the Da operator serves as the marker. This completes the proof. 



D Removing Delta Functions 



In this Appendix, we show how to cancel the delta function factors in (3.28) to obtain the primed Ward 
identities (3.30). To begin with, it is convenient to rewrite the right-hand side of (3.28) as 



limM, 



Qn 



j£o- 



%i ■■■dqi^ \P-yiq) 



1 ^ ^ S^^o 

{Y'<^iq)Oiku . . . , kN))c + i7WT^{Ciq)0{ku ..., kN))c 



^PdQ) 



N 



gn 



+ 



rC^ 



lim M,,„...,„ <j X: [ S^'^ Q^a^ . . . Qj^a^ . ^ ^ 1 Q^a^ . . . Q^a^ J 



— \{0{ki, . . . ,ka + q, . . . ,kN))c 



a=l 



M 



^Tfi^0Wa(^^ 



QT. 



a=l 

N 



dk'i ■■■ dk 



— {0^{ki, ..., ka-l, ka+i, . . . kuhiajaika + q)0'^{kM+l, ■■■, kN))c 



- E ^\u'Hh) ^" 



%3b 



b=M+l 

+ ... , 



dk''. ■■■dk'', 



-{OHh, . . . , 4/)07^,^,,,,^,,...,fc,4,...,^,-^(4/+i, . . . , 4 + g"*, . . . , kN))c 



(D.l) 



where as before the ellipses stand for higher-order tensor terms. In this form, the correlators on both 
sides of the identity involve the same delta function 6'^{P), where P = q + ki -\- . . . k^. 
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First, consider the case n = 0. Expressing (D.l) in terms of primed correlators in this case gives 
= -l_imM,,„HiV<5^^o+f;fc^^j(0(^i,...,%)>'c 

M 

- J2 T''°'"^"(fca)(0^(^l, . . . , 4-1, 4+1, . . . 4f)7*„i„(4)0^(4f+l, . . . , kN))c 
a=l 

N >| 

fe=Af+l J 

-^mM,e,{0{h,...,kN)y,P^^6'{P) + ... (D.2) 

Integrating the last term by parts, and using the fact that P'^6'^{P) = 0, we obtain 

= -M,e, I UN- l)S''o + ^ ^. _|_ j ^o(4, . . . , fe^))'^ 

M 

- 5^ r^°''"^''(4)(0^(4, . . . , 4-1,4+1, . . . kMH,jAka)0^{kM+l, . . . , kN))c 
a=l 

N >| 

- E r^^, >^^(4)(o^(4, . . . , 4.)oL,.,M+....,M.......(^^M+i' • ■ ■ ' ^^^))c + • • ■ 

fe=M+l J 

(n = identity) , (D.3) 

where we have canceled the delta functions from both sides. Note that the limit g — t- has been explicitly 
taken on the right-hand side, hence the right-hand side is now independent of q. This is the desired 
"primed" Ward identity for n = 0, in agreement with (3.30). 

Next we prove that (D.l) holds in the same form for the primed correlation functions for all n > 1: 
q-^O oqe^ ■ ■ ■ dqe„ yP^iq) 3P^(g) J 

^ an 

_ ^ T^^ow.(^j_____(oC(^^^ . . . ^ 4-1,4+1, . . . kMH^jMo^kM+u • . . , kN)yc 

a=l "1 *n 

N ^n -] 

- E r\>^-(4) , . , (0^(4, . . . , ^^m)o7,,^,,,,^,„,.,,,,,,.„.^,^(4.+i, . . . , k^)y. + . . . 

b=M+l "'^h ^'^^n ) 

(n > 1 identities) . (D.4) 
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We prove this by strong induction. (To simplify the notation, we will omit the momentum dependence 
except where necessary, and the g — )• limit will be understood everywhere.) Starting with n = 1, 

\ a \ «i *:o ti / a=l *-i 

^bjb aub \ «M+iiAf+iv,fcb4v«iviiv'c f '^ ^ 

fe=A/+l '^'^^^1 J 

a 



Z^ «6J6 \ iM+lJM+l,---,kb(-b,---iNJN'C ( Qp. ^ ' 



=M+1 



5^'«-"(°>=^"iw5i^*"<^> + - 



,2 \ A^ 



-^^^A E '^^'"^Ir + k^^T;;^ (o);-^r^o->^(o«(fc^i,... ,4-1,4+1,. ..M7w.o^); 



«6J6 /:)l6 \ ^iM+lJM+l,---,kb(-b,---iNJN'C f '^ V'' ^ 

>,=M+1 '^'^^l J 



- Ma,,, I ( (iV - 1)<5^^» + Yl K^] {Oy, - E '^'''''''''{OHki, . . . , 4-1,4+1, . . . 4/)7i„,.0^)c 

I V a ^0 / a=l 

fe=A/+l J ^^1 

where in the last step we have integrated the d'^ /dPi^dPi, once by parts using the fact that Mig^i-, 
is symmetric in its last two indices. Collecting all Qp—6^{P) terms, we see that their coefficient is 
proportional to the n = primed identity (D.3) and thus vanishes. Hence we are left with 



hm M,,„,,-^ ( -l-(y^0(g-)O(4, . . . , ^7v)>'e + -^AmOih, . . . , ^7v))'e 



N / r. . 02 



M „ 

o 



Y r'°*'"^"(4)7^(o^(4, . . . , 4-1, 4+1, • . • 4/)7^...(4)o^(4/+i, • • . , %)); 



a=l *^i 

■ E ^%b''''(f^^)^(^'Ch, . . . , 4.)o7.,,,,.,„...,M...™(^^*^-M' • • • ' ^^^))c + • • ■ 

b=M+l "'^h ) 

(n = 1 identity) , (D.6) 
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which is the desired result for n = 1. 

By strong induction, we next assume that (3.30) holds for all m < n — 1 for some n > 2, and 
show that it must therefore hold for m = n. Expressing (3.28) in terms of primed correlators, clearly the 
desired result corresponds to no derivatives acting on the delta functions. The proof therefore boils down 
to showing that all terms with at least one derivative on a delta function cancel out. Let us first collect 
all ■^pj6^{P) terms, where 1 < j < n — 1. Since M^^^j^ is symmetric in its last n + 1 indices, we get 



Mu,. 



dU^iP) 



n 



Qn~j 



dPi .^ • • • dPe 



Using the fact that 



this reduces to 



n 



Mi, 



j I dqe^ ■ ■ ■ dqt„_^ \ B 



ir'o)', + ^{co)'. 



a 

M 

a=l 

N 

- E r 



n 



gn-j 



+ 



hi 



I n+1 



3Pc 

Qn-j+l 



9kl---dkl_^^ 



{OYc 



n 



Qn-j 



iio ktib 
ibjb 



j J dkl ■ ■ ■ dkl^^ 



=M+1 



n 
3 



{O'^ih, ..., ka-l, 4+1, . . . kMhi^j^O^Yc 



dk^ ■ ■ ■ dk^ ^^ ^M+iJM+i,---,kb^b,---iNJN''^ 



n + l 
j 



n 
3 



n+l 
n — j + 1 



(D.7) 



(D. 



d^6^{P) 



dPf 



n — j + l 



E 



M 



ap,:kr^ i^<^"'"°>^-5«°'' 



^iio 



Qn-j 



ki 



Qn-j + l 



dk'} ■■■dkf 

■t 1 Zfi — j 



\ ^ J^iioiaja 



n 



n — j + l dk*^ 

Qn-j 



dk'^ 



{OYc 



a=l 

N 



j J dkl^ ■■■dk 



n—j+l _ 



E r 



i£o fcb4 



Qn-j 



. ,r . *'^' dk'i ■■■dk^. 

b=M+l h In- J 



{O^O, 



iM+lJM+l<--;kbib,---iNJN''^ 



(D.9) 



Thus the terms with j derivatives on the delta function, with 1 < j < n — 1, are proportional to the 
n — j primed Ward identity, which holds by assumption. It remains to show that terms with at least n 
derivatives on 6^{P) also cancel out. We will see that these are in fact proportional to n = primed 
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identity. Collecting such terms, we obtain 



M 



i^O-^n' 



dPf. ■ ■ ■ dP,, 



^(7'^«o)'. + ^(co);+U.-o + ^,._j^0). 



d 

dk 



M 



E Tr*'°*"'"(0^(^i, • • • , 4-1,4+1, . . . 4/ )7..,.0^)'c 



a=l 

N 



ST^ pj^o ktib/QCQ-y 



iM+lJM+l,---,kb(-b,---iNJN'<^ 



+ - 



1 



b=M+l 
Qn+lj3(p) 



(0>c 



n + 1 ap^o • • • ap^„ ^ 

Integrating the last term once by parts, using the symmetry properties of M, we obtain 



(D.IO) 



M, 



ilO-.^n 



dPf, • • • dPf, 



1 r^"0 / (9 \ 



M 



=1 

N 

ibjb \ iM+iJM+l,-,kb^b,-iNJN'(^ 



a=l 

N 



b=M+l 



(D.ll) 



As advocated, we recognize the n = identity in square brackets, hence terms with at least n derivatives 
on the delta function cancel out. This completes the proof of (D.4). 

E Component Operators 

In this Appendix, we describe a brute force method for constructing the operators used in Sec. 3.4. We 
are looking for operators Pi£Q...i„jmo...m„iQ)i composed from qi and 6ij, which satisfy the four conditions 
listed in Sec. 3.4. 

We start by writing down by writing down all possible terms constructed from qi and 5ij which 
have 2n + 2 indices. First there is a one term which is just a product of 2n + 2 q's, then there are terms 

(2n + 2 \ 
of these since we may 

choose any of the 2 indices for the 6), then there are the terms with two (5's, and so on, until we get to 
the terms which are products of n + 1 (5's. 

Now we take all these terms and write them down with an arbitrary coefficient in front of each. 
Then we impose the conditions listed in Sec. 3.4. Since the conditions are linear, this will generate a set 
of linear equations which the coefficients must satisfy. There will be some subspace of solutions to this 
linear system. Choosing a linearly independent basis of this subspace will yield our component operators. 
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For example, in the simplest case n = 0, we obtain the following 3 operators, 

QioQrno^ij Q£oQj^imo Q.iQmo^C.(,j -* imo^oj ' QiQj^io'mt) ~r ^ij^(.om(, > V^'^j 

where Pij = Sij - qiqj. 

For n = 1 we obtain the follow 3 operators: 

— ^Pimi^loj^limo + ^QiQjPlomiPhmi) ~ ^-^ii^'omi-nimo + ^QloQmoJ'^ijPhmi ~ '^^loQj^imoJ'^lii 



' r^Q-h^mi^imo^lfij r)1hQrrH)J~^imii^^loj r^lh^'mi^ij^lorrnt \ ^Q.hlj^i'mi^lomQ ^^ii^mo Jtjr-»7o"T'i 



mo 
ni 

momi 



' T^Qh^j -'~^i'n'H)-'~^lom.i ~ 'T^QloQrriiJ'^imoJ'^hj ~ ^^/o^mo-Mmi-nijr' + ;^^i9mi-nomo-Mij ~ T^jmi-nomo^'l, 
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-PijPkyrnoPhnn " P^^PlohPj^o + l^PihPhm.Pjmo + l^Pil.Ph^.Pjm, " Pi^,Pl,l,P,m, 

+ ^PUlPlornoPjn.^ + ^PUoPhrnoPjm, + kPijPlohPmom, " Pil.Pl.jPn.,^, " Pu,Pl,jPm,m, ■ (E.4) 

At each n, we will find 3 operators, though their form becomes increasingly complicated as n goes up. 

The operators we find in this manner will not in general be projection operators, because they do 
not satisfy the conditions PaPa = Pa^ab that a set of projectors Pa should satisfy. The subspace of M's 
satisfying the required symmetry conditions has some dimension, d. If the P's were projectors, the sum 
of the dimensions of the images of the P's would equal d, but in our case this sum can be > d, so we have 
an overcomplete set of operators. We can in principle enforce the conditions PaPa = Pa^ab to pare this 
down to a complete set, but since this condition is not linear in the projectors, in practice it is simpler to 
work with the overcomplete set, and we will not be missing anything by doing so. 
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